翻訳と辞書
Words near each other
・ Chirotherium
・ Chirothrix
・ Chiroti
・ Chiroubles
・ Chiroxiphia
・ Chirp
・ Chirp (disambiguation)
・ Chirp (magazine)
・ Chirp (Modern Family)
・ Chirp (scratch)
・ Chirp compression
・ Chirp mass
・ Chirp spectrum
・ Chirp spread spectrum
・ Chirp transmitter
Chirp Z-transform
・ ChiRP-Seq
・ Chirpan
・ Chirpan Peak
・ Chirped mirror
・ Chirped pulse amplification
・ Chirpin'
・ Chirping
・ Chirping cisticola
・ Chirplet transform
・ Chirpy
・ Chirpy Burpy Cheap Sheep
・ Chirpy Chirpy Cheep Cheep
・ Chirpykty
・ Chirpăr


Dictionary Lists
翻訳と辞書 辞書検索 [ 開発暫定版 ]
スポンサード リンク

Chirp Z-transform : ウィキペディア英語版
Chirp Z-transform

The Chirp Z-transform (CZT) is a generalization of the discrete Fourier transform. While the DFT samples the Z plane at uniformly-spaced points along the unit circle, the chirp Z-transform samples along spiral arcs in the Z-plane, corresponding to straight lines in the S plane.〔(A study of the Chirp Z-transform and its applications ) - Shilling, Steve Alan〕〔http://www.mathworks.com/help/signal/ref/czt.html〕 The DFT, real DFT, and zoom DFT can be calculated as special cases of the CZT.
Specifically, the chirp Z transform calculates the Z transform at a finite number of points zk along a logarithmic spiral contour, defined as:〔http://prod.sandia.gov/techlib/access-control.cgi/2005/057084.pdf〕〔
:z_k = A\cdot W^, k=0,1,\dots,M-1
where ''A'' is the complex starting point, ''W'' is the complex ratio between points, and ''M'' is the number of points to calculate.
==Bluestein's algorithm==

Bluestein's algorithm〔(【引用サイトリンク】title=Bluestein's FFT Algorithm )〕 expresses the CZT as a convolution and implements it efficiently using FFT/IFFT.
As the DFT is a special case of the CZT, this allows the efficient calculation of discrete Fourier transform (DFT) of arbitrary sizes, including prime sizes. (The other algorithm for FFTs of prime sizes, Rader's algorithm, also works by rewriting the DFT as a convolution.) It was conceived in 1968 by Leo Bluestein. Bluestein's algorithm can be used to compute more general transforms than the DFT, based on the (unilateral) z-transform (Rabiner ''et al.'', 1969).
Recall that the DFT is defined by the formula
: X_k = \sum_^ x_n e^ nk }
\qquad
k = 0,\dots,N-1.
If we replace the product ''nk'' in the exponent by the identity
:n k = \frac + \frac + \frac
we thus obtain:
: X_k = e^ k^2 } \sum_^ \left( x_n e^ n^2 } \right) e^ (k-n)^2 }
\qquad
k = 0,\dots,N-1.
This summation is precisely a convolution of the two sequences ''a''''n'' and ''b''''n'' defined by:
:a_n = x_n e^ n^2 }
:b_n = e^ n^2 },
with the output of the convolution multiplied by ''N'' phase factors ''b''''k''
*
. That is:
:X_k = b_k^
* \sum_^ a_n b_ \qquad k = 0,\dots,N-1.
This convolution, in turn, can be performed with a pair of FFTs (plus the pre-computed FFT of complex chirp ''b''''n'') via the convolution theorem. The key point is that these FFTs are not of the same length ''N'': such a convolution can be computed exactly from FFTs only by zero-padding it to a length greater than or equal to 2''N''–1. In particular, one can pad to a power of two or some other highly composite size, for which the FFT can be efficiently performed by e.g. the Cooley–Tukey algorithm in O(''N'' log ''N'') time. Thus, Bluestein's algorithm provides an O(''N'' log ''N'') way to compute prime-size DFTs, albeit several times slower than the Cooley–Tukey algorithm for composite sizes.
The use of zero-padding for the convolution in Bluestein's algorithm deserves some additional comment. Suppose we zero-pad to a length ''M'' ≥ 2''N''–1. This means that ''a''''n'' is extended to an array ''A''''n'' of length ''M'', where ''A''''n'' = ''a''''n'' for 0 ≤ ''n'' < ''N'' and ''A''''n'' = 0 otherwise—the usual meaning of "zero-padding". However, because of the ''b''''k''–''n'' term in the convolution, both positive and ''negative'' values of ''n'' are required for ''b''''n'' (noting that ''b''–''n'' = ''b''''n''). The periodic boundaries implied by the DFT of the zero-padded array mean that –''n'' is equivalent to ''M''–''n''. Thus, ''b''''n'' is extended to an array ''B''''n'' of length ''M'', where ''B''0 = ''b''0, ''B''''n'' = ''B''''M''–''n'' = ''b''''n'' for 0 < ''n'' < ''N'', and ''B''''n'' = 0 otherwise. ''A'' and ''B'' are then FFTed, multiplied pointwise, and inverse FFTed to obtain the convolution of ''a'' and ''b'', according to the usual convolution theorem.
Let us also be more precise about what type of convolution is required in Bluestein's algorithm for the DFT. If the sequence ''b''''n'' were periodic in ''n'' with period ''N'', then it would be a cyclic convolution of length ''N'', and the zero-padding would be for computational convenience only. However, this is not generally the case:
:b_ = e^ (n+N)^2 } = b_n e^ (2Nn+N^2) } = (-1)^N b_n .
Therefore, for ''N'' even the convolution is cyclic, but in this case ''N'' is composite and one would normally use a more efficient FFT algorithm such as Cooley–Tukey. For ''N'' odd, however, then ''b''''n'' is antiperiodic and we technically have a negacyclic convolution of length ''N''. Such distinctions disappear when one zero-pads ''a''''n'' to a length of at least 2''N''−1 as described above, however. It is perhaps easiest, therefore, to think of it as a subset of the outputs of a simple linear convolution (i.e. no conceptual "extensions" of the data, periodic or otherwise).

抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)
ウィキペディアで「Chirp Z-transform」の詳細全文を読む



スポンサード リンク
翻訳と辞書 : 翻訳のためのインターネットリソース

Copyright(C) kotoba.ne.jp 1997-2016. All Rights Reserved.